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Abstract. Intertwining relations for A^-particle Calogero-like models with internal degrees 
of freedom are investigated. Starting from the well known Dunkl-Polychronakos operators, 
' we construct new kind of local (without exchange operation) differential operators. These 

operators intertwine the matrix Hamiltonians corresponding to irreducible representations 
Q>^ . of the permutation group Sjsr- In particular cases, this method allows to construct a new 

I class of exactly solvable Dirac-like equations and a new class of matrix models with shape 

I invariance. The connection with approach of multidimensional supersymmetric quantum 

' mechanics is established. 

p i' 

(D ! 1. Introduction 

> 

^ . The exactly solvable quantum A^-body problems have provided useful tools to 

^ I investigate both formal algebraic and analytic properties with applications to different 

branches of Physics. The most intensively studied models are the Calogero model 
(many-body extension of the one- dimensional singular harmomic oscillator) [||], 
and its various generalizations, so called Calogero-like models. The latter ones have 
either scalar 0-0 or matrix (with internal degrees of freedom) 10- nature. The 
Calogero-like models have been widely developed incorporating many-body forces 
[ p!0[| , different root systems |TT[,|T^ and mult i- dimensions |J13[. The supersymmetric 



extensions of Calogero-like models ||I5-|[T^ also seem to be very promising. 



In the papers [0, [gg, §, |2T[, |22[, § different types of Dunkl operators 
for the investigation of Calogero-like models were usedi. These operators intertwine 



°We call the one-particle operators constructed in the Dunkl-Polychronakos operators 

(DP) to distinguish them from the genuine Dimkl operators set forth in the papers p^ , 
1^, which are slightly different. 
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Calogero-like Hamiltonians and therefore allow one to construct the integrals of mo- 
tion and the eigenfunctions (if they exist) for these models. 

The characteristic trait of DP operators and of the corresponding integrals of 
motion derived from them is that they involve the coordinate exchange operators, 
and thus are nonlocali. These techniques are briefly outlined in Subsection 2.1 of the 
paper. 

On the other hand, the multidimensional supersymmetric quantum mechanics 
(SUSY QM) [Q, applied to the Calogero-like models |]l3[-|[T^,[Q provides one with 
another set of the intertwining relations, where the matrix Calogero-like Hamilto- 
nians of a specific type 



18 



are intertwined by the the supercharge operators. In 
that approach, both the Hamiltonians and the supercharge operators are local. The 
intertwining relations are the most important part of the SUSY QM algebra, which 



25 , 26 



is clear from a number of generalizations of the standard SUSY QM, e.g. 

In the rest of Section 2 these two approaches will be unified: from the DP op- 
erators set forth in pO[, we construct new local operators of the first order in 
derivatives, that play the role of intertwining operators between the matrix Calogero- 
like Hamiltonians. For the Calogero-Sutherland model [Q] this leads to a new class of 
exactly solvable Dirac-like (matrix and of the first order in derivatives) Hamiltonians. 
For the Calogero model with oscillator terms the new intertwining relations give a 
new implementation of the shape invariance condition |2^,[^, ||T5[|,|p6||. 

In Section 3 we consider the particular case of three-particle Calogero-like models, 
which is the simplest nontrivial realization of the method introduced in Section 2. 
For the models without oscillator terms (OT) the above mentioned Dirac-like Hamil- 
tonians coincide with the conventional Dirac Hamiltonians for a massless particle in 
magnetic field. In that case, the 2x2 matrix Calogero-like Hamiltonians can be 
interpreted as the Pauli Hamiltonians for the same system. 

In Section 4 the class of the Calogero-like matrix Hamiltonians described in 
is considered. For such models the intertwining relations derived in Section 2 are 
reduced to the well-known SUSY QM relations [^. However, there is a wide class of 
models for which the intertwining relations introduced in Section 2 are not reduced to 
any previously known ones. Clearly, the SUSY QM is valid not only for the Calogero- 
like models, but for many other multidimenional and multiparticle ones [^]. The 
question as to how far the generalization of SUSY QM constructed below can be 
extended^ to non Calogero-like models, deserves further attention. 

A possible way of extension of the formalism presented below is to consider the 
generalizations of the Calogero-like models incorporating many-body forces [|l^ and 



^By the words "nonlocal" and "local" we mean "containing exchange operators" and "containing 
no exchange operators" . 

'^One example of such extention (for the three-particle case) will be given in the first footnote on 
the page 11. 
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different root systems |T^, for which the Dunkl operators also exist [Tl[],|]T2|, and 
therefore one can construct from them local intertwining operators analogous to those 
of the present paper. 

2. Intertwining operators of first order in derivatives. 

2.1. Dunkl-Polychronakos operators [^Df| ,[|6|1 for Calogero- 
like models. 



Let us consider a one-dimensional quantum system of N particles with coordinates 
Xi- Let Mij be the operator that exchanges the coordinates of the i-th and j-th 
particleJi. The DP operators are defined ||2^,[§[ as: 



TTi = -idi + i ^ VijMij = Tcj Vij = V{xi - xj); (1) 
^^ = ^■, V{x) = V{-x) = V*{x). 

OXi 

The operators vTj are one-particle ones, i.e. 

Mi.Tii = TTjMij; [Mi,, TTk] =0, I, J. (2) 

Their commutators can be written as: 

[^i'^i] = Yl Vijk[Mijk - Mjik], where 
Vijk = VijVjk + VjkVH + Vk^Vij, (3) 
and Mijk are the operators of cyclic permutations in three indices: 

Mijk = MijMjk = Mjki = Mkij = M^^. 

In the cases both of the Calogero-Sutherland (CS) models 

V{x) = /cot a; (trigonometric or TCS model), (4) 
V{x) = /coth (hyperbolic), 

and of the delta-function model 

V{x) = I signx, (5) 



"^The small latin indices range from 1 to iV everywhere . 



4 



the function Vijk = so that 



(6) 



The Hamiltonians for these models are 
H 



(7) 



where A = J2ididi and V-j = V'{xi — Xj). It is known pO|,||§l that in this case 

KH] = 0. 
In the case of the Calogero model, 



V{x) = l/x] 



ijk 



0; 



0, 



(9) 



and the equations (0),(|^) remain valid. What is more, the DP operators themselves 
mutually commute. However, this model doesn't have a discrete spectrum. 

The Calogero model is usually considered in a harmonic confining potential (we 
abbreviate it as CO: Calogero-Oscillator). For this model, the following operators 
should be introduced ^M, (see also fllj, [pl|],[|22|,|]): 



aj = TTi ± lujXi; 
The Hamiltonian can be written as: 

Hco = E 



a+a- + IcoJ:M., = -A + c^E^? + E + 

[Xt Xj ) 



(10) 



It has been proven that the operators Hco, o,f form the oscillator algebra: 

[Hco,af] = ±2uaf. (12) 



'^The following procedure is applicable to arbitrary set of operators tt^, provided they satisfy 
(|^),(^ and the Hamiltonian is given by the second equality of (|^). 
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2.2. The local form of the Hamiltonians. 



Let us consider an irreducible representation A of the permutation group Sn 
realized on real vector functions /^(xi, ...,X]\f); a = 1, ...,dimA by the matrices T^-: 

M,,/, = iT,^)paU, (13) 

where {T^)pa = {T^)ap are the (constant) matrix elements of the permutation oper- 
ator Mij in the representation A. Below we will assume summation over the repeated 
indices, unless specified otherwise. We will also use the fact that are real 
symmetric orthogonal matrices. 

It will be useful to introduce the vector notations: 

f = e„/„, (14) 

where the constant vectors e^, [a = l,...,dimy4) form a basis in the space of the 
representation A. Then it is also helpful to define the operator in the vector form: 

Tge„ = e^ie^yTf^e^ ^ {T.^Uep. (15) 

Multiplying ( P^ onto and using (|T^, we get: 

M,,f = Tgf, (16) 

where M^j act only on the arguments of fa and Tf- only on the vectors e^. 

All the Hamiltonians H from the previous Subsection can be written a.s H = 
Hscai + where Hscai are scalar operators containing no exchange operator 

terms: 

H,,,i = -A + Y.V^^ (17) 
for the models without OT (|)-(|), and 

for the CO model (0). The Hamiltonians H act on the functions0 from the repre- 
sentation A as 

Hfa = Hpjp] H^^ = HscalSal3 + ^iji'^ij)f>» = ^^0, (19) 



^Note that the functions f are not necessarily eigenfunctions of H . We do not discuss in the 
present paper the symmetry properties of the eigenfunctions of the Calogero-hke models. 
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or, in the vector form, 

m = H^f . 

Note that if f satisfies (|TB|) then Hf satisfies it too, since Eq. (^) is equivalent to 
the condition: T^Mj^f = f (no summation over The latter condition is satisfied 
for Hf, because [H, TfjMij] = 0. 



2.3. The intertwining operators in the local form. 



The matrix Hamiltonians (19) do not contain exchange operators Mij explicitly. 
Our aim now is to get rid of the Mij in the DP operators (|ip,(|TD|) too, and rewrite 
the Eqs. (@),(|l3) in terms of local operators only. 

Let us study the action of the DP operators on the symmetric functions satisfying 
(p!3D,(p!6D. The expression Hifa no longer satisfies (|13D even when fa satisfies it. 



Instead, iiifa transforms under the action of M^j as an object from the direct product 
of representations for vTj and fa- Of course, the DP operators transform under Sn 
in accordance with (|]). However, the Tr, belong to a reducible representation of Sn 
because tti + ... + tcn realizes the absolutely symmetric representation. Therefore, it 



is helpful to go to the well-known Jacobi coordinates ||29|| : 



= I ^ (xi + ... + xg-ea:g+i); 1 < ^ < - 1 (20) 

+ 1) 



1 ^ 



ori Uk = RkmXmi where the orthogonal matrix R is determined by ([20|). The deriva- 
tives are connected by the same matrix: d/dyk = Rkmd/dxm, because R is an or- 
thogonal matrix. Similarly, we can write the DP operators in the Jacobi variables: 

Pi = I / (7ri + ... + 7rg-e7rg+i); 1 < ^ < - 1 

ve(e + 1) 
1 ^ I 



fThe indices of the Jacobi variables denoted by Greek letters range from 1 to N-1; those denoted 
by Latin letters range from 1 to N 
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or pk = RkmT^m- The operators now transform under Sn as the irreducible repre- 
sentation r with the Young tableau (A^ — 1,1). Similarly to (pTsD this fact can be 
written as: 



Mr,p^ = {Tf^x^px. (21) 



This is a property of the Jacobi variables (see e.g. [0). 

The object p^fa transforms under the action of Sn as the interior product V x A 
of the representations A and F, or, in more detail, in accordance with the formulae 
(0),(^, as 

Mijp^fa = {Tfj) x^{Tf^) paPxf (3- 

As outlined in the book ||30|(chapter 7,§13), the interior product V x A contains only 
the irreducible representations of S'at, whose Young tableaux differ from the tableau 
for A by no more than the position of one cell (but not necessarily all of them). For 
example, for the absolutely symmetric representation of with the Young tableau 
(A^), obviously, 

F X (iV) = F, 

and the result does not contain (A^). 

Let B be some irreducible representation that appears in F x A. Then we can 
extract its contribution to F x A with the help of the Clebsch- Gordan coefficients 
{Vi,Aa\Ba) = (^a|a): 



= {^a\a)p^fa = D^afa] Daa = ({a|cr)pg. (22) 

The resulting function g^^ satisfies the analog of (|13D for the representation B: 

M,,g^ = iT,^)sags- (23) 

This can be checked directly by substituting 
expression!: 

{T[^)xdT.j)p^{^a\a) = (7;f ),.(A/?|5), (24) 
and the fact that T,.,- are hermitean matrices. 



''The standard notation for a Young diagram containing Xi cells in the i-th line is (Ai, . . . , A„); 
if the diagram contains m identical lines with fi cells, it is denoted by (. . . , /x™, . . .). 

'The expression (|4|) (see |3^, formula (5.114)) is actually a necessary and sufficient condition for 
(C^, Aa\Ba) to be a Clebsch- Gordan coefficient for arbitrary representations A, B and C & Ay. B. 



22) into (B3) and making use of the 
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On the functions that satisfy (|1^) the operator D^a acts as: 

Daa fa — /q , 

= m^)ipt)po. = {m^)R^k{7rt)pa, (25) 

(vr^)/3a = -idkSf3a+iJ2^km{TkJ(3a- (26) 

For the models without OT [H, vTj] = 0, so, [H, D^a] = 0. Hence, for any 

HD^afa = D^aHfa- (27) 

For all fa which satisfy (|I^), H^^fp satisfies (|I^) (see the end of the previous Sub- 
section), and D„afa satisfies (pS]). Using these symmetry properties, we can expand 
the sides of the equation (^) as 



(l.h.s.) = Hf^Dsafa = HlDfj, 



on 



(r.h.s.) = D^aH^Jp = D^Hlfp. 

Taking into account that H"^, are symmetric matrices in the internal indices, we 
can see that on the functions satisfying (|13D 

H^sK = DtHU (28) 



Similarly to the above, for the CO model from the Eqs. (|1^) it follows that 
HD^^ = = D^^{H ± u). Following the same route that led us from (pTl) to (^8]) , 
we can conclude that on the functions satisfying (p!3|), 

Hi^Dt^ = D^^ {H^p ± 2uj5ap) , (29) 

where 

Dta = m^)RiAi4)pa±^uJX,6ap]■ (30) 



Note that all terms in the Eqs. (PS]),(EU|) are local: they contain no exchange operators 
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2.4. The operatorial nature of the intertwining relations. 



The Eqs. (P8D ,(|29[) are not yet operatorial intertwining relations such as, for 



example, the SUSY QM ones (see [^,^T|,[^), because the former are valid only 
on the functions that satisfy the symmetry condition (0). On the functions outside 
that class the Eqs. (pSf) , p9|) , generally speaking, may no longer be satisfied. 

However, we can prove, that they are satisfied on all functions and thus are oper- 
atorial intertwining relations, by making use of the following 

Theorem 1: Let A be some representation of Sn- Let Lap he some linear dif- 
ferential operator of finite order with the coefficients being rational matrix functions 
of the variables Xi, or sinxj,cosXi, or sha;j,chxj (but not of any two of them simul- 
taneously), singular a^ U = {x|3i, j : i ^ f Xi = xj} at most. The coefficients are 
matrices of dimension dimA x dimA. Then, if 

Lapfp = 

for all fp satisfying ([7^, then L = as an operator. 

The proof of this Theorem can be found in the Appendix 1. Using this Theorem 
1 for the difference of the left and right parts of the Eqs. we can conclude 



that the latter are =atMed operatoriall^. 

In particular, it means that that when the initial representation A coincides with 
the resulting representation B, the operators are integrals of motion for the 
(trigonometric and hyperbolic) matrix CS models. Therefore, each CS matrix model 
corresponding to a representation A such that A & T x A has a local integral of 
motion D"^ of the first order in derivatives. An example of model from this class will 
be given in the Section 3. Note that e.g. the models with A = (N) or A = (1^) lie 
outside this class. 

With periodic boundary condition Jon a unit circle S, the Hamiltonians H"^ for the 
TCS system are exactly solvable (see e.g. 0) and have discrete spectrum and finite 
dimensional degeneracy of levels. The fact that (^), ( P^ ) are satisfied operatorially 
when B = A allows us to find also the spectrum and the eigenfunctions of the D"^ , 
i.e. the normalizable functions /^(x),^ = 1, ...,dimy4: 

D%fp = efa- 



■'By X without an index we mean the vector (xi, ...,xn) with N elements. 

'^The Theorem 1 cannot be apphed to the model with V{x) = Zsignx, because the Eq. (|2^ ) may 
then contain delta- functions-like singularities at U . 

'it means that Xi £ S = [0,2tt], or x e 5^, where S'^ is a torus; {{xi,...,Xi + 2tt,...,xn) — 
i{xi,...,x.„ ...,xn)- 
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The operators D"^ for the TCS system may be considered here as Dirac-hke Hamil- 
tonians of first order in derivatives: 

/^From the commutation relations (p8|) it follows that for a given A, the operators H"^ 
and D"^ can be diagonalized simultaneouslyH. In more detail, if f^^"'* are degenerate 
eigenstates of H"^ with energy En- 

n^t^ = Ejt\ (31) 

and b = 1,..., J (J is the degree of degeneracy), then after acting by on both 
sides of the equality (|3T|) we see that D^4^"^ satisfies (0) too. HenceB, there exists a 
constant J x J matrixO _F("): 

Because H"^ is hermitean, we can choose 4^"^ that constitute a basis of periodic 
vector functions with dim A components on . Because D"^ is hermitean, one can 
check that F'-"-' is also hermitean, and it can be diagonalized by a unitary rotation: 

'^ab ^bc '-^dc — "ac-, 



and the functions 



>^(x) = f/rfr^(x) 



are eigenfunctions of the operator D"^ with energies e^"-*. Because L^^J^"* are unitary 
matrices, the functions g^"-* constitute a basis of the periodic vector functions on . 
Hence, they form a full set of eigenfunctions of D"^. 

For the CO model the case of S = A is interesting too. The Eqs. ( pOj) will then 
take the form: 

HtDf^ = D^^ {H^p ±2u^5^p), (32) 
where D"^^ are still defined by (|30|) , and will be satisfied operatorially. They are the 



relations of the oscillator- like algebra (similar to ([T^). In the language of SUSY QM 



™See [^; compare also with the SUSY QM intertwining relations which were used to find a part 
of spectrum of Hamiltonians in one jsj] or two |Q dimensions. 

"We assume for simplicity that the action of D"* does not destroy the normalizablilty of f^'""* . 
°Note that the matrix elements F^'f' are NOT dimyl x dim A matrices, but just scalar constants. 

(n) 

In other words, they do not affect the vector structure of . 
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it means that each CO matrix model for a representation A ^ T x A obeys shape 
invariance (SI) in A^ — 1 dimensions (because the center of mass motion is decoupled). 
Let us stress that this is the first example of SI in several dimensions realised by local 
operators of the first order in derivatives (cf. the attempts in [0). The nonlocal SI 
of the Calogero model (see the Eq. (|12D) was described in |2D,[§, lll[, [0, [0, §; 
a model with a two-dimensional SI of the second order in derivatives was proposed in 



3. Examples with N = 3. 



Assume that = 3, and both representations A and B are equal to F = (2, 1). 
Then the functions that satisfy ([T6|) are: 

/211 
, /121 . 

where the index a = (211), (121) enumerates the partitions of the Young tableau 
(2, 1) (see or the end of chapter 1): 



P+/121 



P /211 — 0; M12/211 — /211; 

P /121 = 0; M12/121 = — /121; 

1 + M12M13 + M13M12 ± (M12 + M23 + M13) 



are the projectors onto the symmetric/antisymmetric representations of 5*3. The 
matrices of the representation F for = 3 have the form: 



■12 



(^3, J-23 = -^^1 - 



(Ti 

2 



2^3, 



(33) 



where cTj are the Pauli matrices. The Hamiltonians for the Calogero-like systems 
without OT for the representation F have the form (p!9|): 



Hr = -A 



E 



y.2 + V'T^- 



(34) 



As to the DP operators in the Jacobi variables, one can check that pi has the 
same symmetry as /121, and p2 - the same as /211. The Clebsch- Gordan coefficients 
for F X F — >■ F are written in the end of chapter 7 of the book Plugging them and 
the T,^^- from (|33| ) into the definition of D'" (expressions (l25|),(^)), we can conclude 
after some algebra that for the Calogero-like models without OT 



1 

V2 



d d 

-ia^- hiai- V2(J2{Vi2 + V23 + V31) 

dy2 oyi 



(35) 
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where are the Jacobi variables (|20|). The Eq. ( 
(|35D is satisfied operatorially, i.e. [H'", D'"] = 0. 
iV = 3 can be proven by direct calculation: 



for from (^) and D'" from 
In fact, a stronger statement for 



fD 



r\2 



dyi 



+ C, 



(36) 



where C is a real constant. What is moreO, it is true operatorially even for the case 
V{x) = I sign X. In that case we were unable to prove (^) for arbitrary representation 
A, but for = 3 and A = T it turns out to be true. 

Eq. (^) signifies that the operator D'" realizes a sorti of a "square root" of H'". 
This, in particular, means that the spectrum and eigenfunctions of itself can 
be found easier than in the general case described in Section 2 of this paper, if the 
spectrum and eigenfunctions of are known. 

For the TCS model the operator (RH) has the form: 



1 

V2 



_d_ 

dy2 



d ( 
— -/V2a2fcot(v^l/i: 



cotf — 



+ cot| 



■y\ - 




The eigenfunctions of this operator are 2-component column functions g(?/i, y^)'- 

D^g(l/i,Z/2) = eg(yi,?/2). (37) 

In this case it follows from (|36|) that g is also an eigenfunction of with energy 
E = 2{^-C). 

Let us prove now that all the eigenfunctions of D'" can be obtained from the ones 
of H^. Let f (x) be an eigenfunction of H'" with energy E: H^f = Ei and zero total 
momentum. Then the following alternative should be considered: 

a) f itself is already an eigenfunction of D^, i.e. it satisfies (^). Then the 
corresponding eigenvalue is: e 



±^E/2 + C. 

b) D'^f = u and f(x) are linearly independent. Then one can check that D^u = 
(C + i?/2)f, and u(x) is also an eigenfunction of with energy E. Thus, the 



PQne could even replace the term V\i + V23 + V31 in ( p5| ) by an arbitrary function w(yi, y-i). Then 
the square of D'" would remain the sum of the Laplacian and a momentum independent 2x2 matrix 
potential. However, we do not know any cases when this sum is an exactly solvable Hamiltonian, 
except for those given in this text. 

''The operator (35) can be viewed as a Dirac operator for a massless fermion in three dimensions 
(2/1,2/21^3) in the magnetic field that does not depend on 2/3 and is orthogonal to the axis 2/3- The 
component of the fermion's momentum along the axis 2/3 should be zero. The Hamiltonian (|3^) is 
then the Pauli Hamiltonian for the same system [p6[,pl[. 
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eigenfunctions of H'", that are not the ones of D'", form pairs in which D'" transforms 
each member into another. If f is normahzable then u is too, because 

< u|u >=< ilCD^fli >= {E/2 + C) < f |f >< +00. 

/^From each such pair one can construct two eigenfunctions of D'": 

f ± (C + E/2)-i/2^ 

with energiell e = ±JC + E/2. 



Thus, when taken in the above form, the sets of eigenfunctions of H + ^ and 



concide. 

For the CO model the Hamiltonian 
form: 



for the representation F will have the 



\Xi Xj) 



where Tjj are defined in (|33D . The intertwining operators (0) can be rewritten as: 



D 



r± 



1 

V2 



-20-3 



_d_ 

dy2 



-V2(T2( — + 



+ 



d_ 

dyi 



V3, 



The operatorial relations (|32D of the oscillator-like algebra 

correspond to the SI of the matrix 2x2 Hamiltonian in two dimensions (because 
the center of mass motion is decoupled). 



4. Connection with the ordinary SUSY QM. 



In this Section we will restrict ourselves to the class of representations with Young 
tableaux of the form 

A = {N-n,ry, n = l,...,N. (38) 

' The situation C + E/2 < is impossible because otherwise the Eq. ( ^ ) would remain valid, 
but with imaginary e. The operator D'" is hermitean, so C + E/2 > 0. 
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It was proven in |]T8|] that for this class of representations we can choose a basis (see 
(|l4D)with the help of the fermionic creation/annihilation operators i = 1...N 

= 0, {^+,^+} = 0, {^„^+} = 6,,; (39) 
ipi\0> = 0; i,j = l...N; < 0|0 >= 1. 



It is useful to introduce also the fermionic analogues 0^ of the Jacobi variables (|2 
(see 0): 

where Rkm are defined in (pOf). The fermionic Jacobi variables obey anticommutation 
relations similar to (IST 



= 0, {0+,0+} = O, {0,,0+} = 4r„. (40) 
0fc|O> = 0; k,m = l...N. 

Now we can define the basisS e,^: 

= 0+ ...0+jO >= > = 1, iV - 1, (41) 

where a = {ai...an) is a multiindex with values in the fermionic number space, and 



GaY^a = < a„...ai|ai...a„ >= 1 



(e^)"^ = {\ai...an >y =< O|0a„---0ai an---ai| (42) 

(no sumation over a is implied). 

Because of (^OD, it is sufficient to include into the basis only the vectors with, 
say, «! < ... < an, and the summation over a will be done over such vectors only. 
It was also proven in [T^, that for the basis (^Tj) the operator Tfj in the vector 



form ( P^ can be realized as 

iT,^)a(sei3 = Tf^Ga = T^|ai...an >= K,j\ai...an >, (43) 



where 



K^, = ijf^i + ^ti^i - ^ti^^ - ^ti^j + 1 = 1 - (^+ - - V^,) = 

= 4, = {k,^)\ (44) 



"The fermionic operators ^^r, (/)^ do not enter into e^, because they correspond to the center of 
mass degree of freedom which is decoupled. 
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It follows from (^) that all the Hamiltonians H"^ (^) with A from the class (|38|) 



take the same form in the basis (|4l|) : 



H = Hseai + J2v:,K.r 



(45) 



Let us consider the intertwining relations ( P8[ ) for the Calogero-like models without 
OT for A from the class From (|5]),(|l3) it follows that the Hamiltonians H^, 
in will have the form: 



H 



(46) 



given m 



One may notice that this Hamiltonian is a particular case of the Superhamiltonian 
]TB|, up to the sign of Vij and an additive scalar constant. 

) and can differ by 



Because the Young tableaux for B and A belong to the class ( 
no more than the position of one cell (see ^0|) chapter 7,§13), B can either coincide 
with A or have the form (A^ — n =]= 1, 1"^^). 

Let us consider the case B = {N — n — 1, 1"^^), for which we can realize (cf. ||ll 
the Clebsch-Gordan coefficients in the operators D"^ as: 



(47) 



One can check that the Clebsch-Gordan coefficients, defined by (P7|), satisfy (^), and 
therefore they correctly connect the representations A = [N — n, 1"), F = (A^ ~ 1? 1) 
and B = {N — n — 1, 1"+^). These Clebsch-Gordan coefficients may differ from the 
standard ones (see |^) by an inessential overall factor. 

Now we can express the intertwining operators D"^ (^5]) in terms of the fermionic 
operators defined above: 



D^e^ = e,(e,)'^D^e^ = e^D 



A 

al3 



\ai...an+i >< an+i...ai\C,ai...an > R^k 

= R^k 



-^^k5pa + ^EyUTkt) 

mj^k 



m^k 



R^k4>^ 



idk\Pl...Pn > +iYl Vkm(T^m)l3o.\(^l---(^n > 



R^k4>^ 



m=^k 

-idk + i VkmKkm 
ra^k 



|A.../3n > . 
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So, we can conclude that 



-idk + i XI ykmKkm 



(48) 



The operator D"^ turns out not to depend on the specific choice of A from the class 
(|38D , i.e. on the fermionic number. It augments the fermionic number by 1; that is 
related to the fact that D"^ changes the position of one cell in the Young tableau for 
A (see ig). 

One can notice that 



1 



(49) 



and after some algebra one can deduce from 

. ^ d 



dy 



N 



+ rk 



(50) 



To further simplify the form of the intertwining operator D"^, one can use the following 
Theorem 2: For all Vkm with Vkm = —Vmk the following equation is satisfied: 

i't VkmKkm = i^tVkm, 
mj^k m^k 

where Kkm is the fermionic permutation operator ([^/ 

The proof of the Theorem 2 can be found in the Appendix 2. 
Now one can rewrite the Eq. (^) as: 



—iq 



dk — 'Y Vkm 

m^k 



(51) 



(52) 



Therefore, the Eq. (pSf ) with the Hamiltonian (^61) takes the form: 

[i/,g+] = 0, 

and its hermitean conjugation! is: 

{H,q-] = Q- q~^{q^)l (53) 

*The Eq. (|5^) could also be obtained in another way: we could consider the Eq. ( p8| ) with 
A = {N — n, 1") but with B = {N — n + 1, 1"~^). Using formulae similar to (^7|)-(|5l|) one can check 
that then = iq^ = i{q+y. 



17 



The operators (^) coincide with the supercharge operatorsS for the Calo- 
gero-hke models given in ||18[, except for the sign of V^j, which is determined by the 
sign of the constant I in Therefore, if we replace I by — / in the SUSY QM 

relations of [0, they can be rewritten as: 



{q-^q+] = H + ^ + C- {q^)' = iq 



0. 



(54) 



The term -^p- in (Q) is unimportant because it commutes with and H. 



The commutation relations (|52| ) , (|53|) may be considered as the SUSY QM com- 
mutation relations, corresponding to the algebra ( |5^ ) with the supercharges q^ and 
the Superhamiltonian H + -^p- + C. 

We can conclude that for the models without OT the intertwining relations (|28|) 
for the representations A, B from the class ( pHf ) turn into the relations of SUSY QM 
, ||18|| . For other A and ^ ( |28D can be considered as a generalization of the SUSY 



, - For other A and B 
QM intertwining relational. 

Now let us turn to the CO model and the intertwining relations (p9|) with A from 
the class (^). From (^),(|18]) it follows that the Hamiltonians H"^,H^ in 



are: 



l{l - K, 



+ Nuj. 



[Xi 



(55) 



This Hamiltonian, analogously to the previous case (0), has the same form for all 
representations A from the class (^). However, the Hamiltonian (^) differs slightly 
from the corresponding Calogero Hamiltonian given in , as will be explained below. 

The intertwining operators D"^^ can be treated similarly to the case without OT, 
the only difference being that one should write di =F uxi instead of di everywhere. In 
particular, for the case with B = {N — n — 1, l""*"^) the definition ( pO|) leads to the 
following analog of the formula (^8]) for D"^^ (the same for all A from the class (p8|)): 



^ km 



Rik<P'l -idk ± ii^Xk + i^{xk- Xm) ^K, 

m^k 



(56) 



"The operator from ( pT|) differs from the standard supercharge operator for the TCS model 
jl7[| by the term ~4'n^^ that cancels the dependence of the supercharge on the center of mass 

coordinates yN, ipN- 

'"The intertwining relations are the most important part of the SUSY QM algebra, which is clear 
from a number of generalizations of the standard SUSY QM: e.g., [E5|,pffl. 
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Making use of the Eq. ( ^9[ ) and of the Theorem 2, we can obtain that, similarly 

to my. 



Wkm 



= -<Pn 



d 



N 



dk 



J2 ^km 



(57) 



W{Xk - Xr. 



W{x) = ^x + -- 
JM X 



A. 

dy 



Wkm = W{Xk-Xm); 



'N7r- + ^k dk - 



'N 



Wix) = X + -. 

N X 



Taking into account the formulae (|^ , , we can rewrite the Eq. 
and A from the class fBSl) as: 



and its hermitean conjugation0: 



[H,q- 



—2ujq 
d 



dvN 



[H,q 
-dk - 



-2ujq'^ 



-] = 2ujq 

m^k 



for 

(58) 

(59) 
(60) 



(g+)t = 



d 



N 



dy 



+ 'ipk 



N 



-9k - J2 ^km 

m^k 



The operators q from (^^ , (|60D are similar to the supercharge operatoril q given 



18[| for the Calogero model with OT, up to a redefinition of constants. Therefore, 

(61) 



m 

we can construct the following SUSY algebra: 

{q-,q+] = h- {q^f = {q-f 
with the Superhamiltonian h: 



0, 



h = H + 2u^^^k - 



N, 



H 



N 



dy 



(62) 



N 



'"The Eqs. ( p9D could also be obtained in another way: we could consider the Eq. (M) with 
A = {N-n, 1") but with B ^ {N Using formulae similar to (|^-(|l]), (||), (f^) one 
can check that then D"^+ = iq^ , where is defined in (|60|), and D'^^ = iq^ . 

^The operator q^ from ( ^7| ) differs fr'om the standard supercharge operator for the TCS model 
by the term ~(t>N'^^ that cancels the dependence of the supercharge on the center of mass 
coordinates hn and (j^jq. 
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where H is defined in (|55|), and C is a scalar constant. The term ifjv is unimportant 
because it commutes with and H. The operators form an algebra similar to 
(|6TD but the sign of u in the Superhamiltonian (|62D should be different (cf. ||37|1): 



0; 



h = H - 2uji)tiJk - H, 



N, 



H 



92 



N 



+ UJ yj^ - 2uj(f)%(f)N + C. 



/^From the SUSY algebrae (0),(^) one can deduce the commutation relations that 



can be shown to be equivalent to (|58D ,(|59D: 



[h,q' 



0: 



[h,q 



;?±1 



0. 



The Eqs. (^),(^) in the case of A,B from the class (pq ) are reduced to the ordi- 
nary multidimensional SUSY QM for the Calogero-like models [0,||ll|- However, 
for A or B outside that class the Eqs. 



(I28|),p9|) describe a generalization of the SUSY 
QM intertwining relations that has not been known before. Clearly, the SUSY QM is 
valid not only for the Calogero- like models, but for many others [^. The question 
as to how far the generalization of SUSY QM constructed above can be extended to 
other, non Calogero- like models, deserves further attention. 
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Appendix 1. 

In this Appendix, we prove the following 

Theorem 1: Let A be some representation of S^. Let L^p he some linear dif- 
ferential operator of finite order with the coefficients being rational matrix functions 
of the variables Xi, or sinxj,cosXi, or shxj,chxj (but not of any two of them simul- 
taneously), singular at U = {x|3i,j : i ^ j, Xi = Xj} at most. The coefficients are 
matrices of dimension dimA x dimA. Then, if 

Lapfp = 

for all fi3 satisfying j^lBj), then L = as an operator. 

Proof: Consider the principal Veyl chamber: {x : Xi < ... < Xn}- Every function 
defined on this chamber can be continued onto the rest of gfj'i™^ by using ([T^). The 
result will obviously satisfy ([TsD , so it is annihilated by L. Hence, L annihilates all 
functions on the principal Veyl chamber. The same can be stated about every other 
Veyl chamber: {x }. Hence, L annihilates all functions on 3?'^'™'^ \ U . 

From the fact that the coefficients of L are rational functions of Xj, or sinxj, cos Xj, 
or shxj, chxj, it then follows that they are zero identically. 



Appendix 2. 



In this Appendix, we prove the following 
Theorem 2: For all Vkm such that Vkm = —Vmk, 

where Kkm is the fermionic permutation operator defined in 

Proof: taking into account the definition (|4^), we can check that 



(63) 



(64) 



In the Eq. (|6^) no summation over either index is implied. Substituting (|6^) into the 
left side of (IBBI) we see that 



mj^k mj^k 



E 

m^k '" 



m^k 
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